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£S) Abstract 

O I We show that the optical structure of the helical phase of a chiral nematic is naturally 

associated with the Bianchi VIIo group manifold, of which we give a full account. The 
Joets-Ribotta metric governing propagation of the extraordinary rays is invariant under 
the simply transitive action of the universal cover E{2) of the three dimensional Euclidean 
group of two dimensions. Thus extraordinary light rays are geodesies of a left-invariant 
metric on this Bianchi type VIIo group. We are able to solve by separation of variables 
both the wave equation and the Hamilton-Jacobi equation for this metric. The former 
reduces to Mathieu's equation and the later to the quadrantal pendulum equation. We 
£h discuss Maxwell's equations for uniaxial optical materials where the configuration is invariant 

i under a group action and develop a formalism to take advantage of these symmetries. The 

material is not assumed to be impedance matched, thus going beyond the usual scope of 
£Nj transformation optics. We show that for a chiral nematic in its helical phase Maxwell's 

^ equations reduce to a generalised Mathieu equation. Our results may also be relevant to 

helical phases of some magnetic materials and to light propagation in certain cosmological 
models. 

\Q 1 Introduction 

O 

Recent years have seen a growing interest in the application of the geometrical ideas originally 
developed for studying Einstein's theory of General Relativity to other areas of physics, such 
as condensed matter physics. The motivation is both the theoretical aim of developing the 
ITj mathematical tools capable of dealing with as wide a range of physical problems as possible, 

and the desire to construct laboratory analogues of the exotic conditions which general relativity 
allows, but which are likely ever to remain inaccessible to direct experimental investigation. This 
in turn may provide a stimulus for further laboratory investigations, possibly resulting in the 
discovery of new physical effects. 

In the present article we shall pursue this direction by demonstrating how the mathematical 
formalism of Lie groups, which is of widespread use in General Relativity and High Energy 
Physics [1, 2], can be harnessed to study optical properties of symmetrical phases of matter. 
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We develop a formalism allowing the symmetries of an electromagnetic medium to be directly 
exploited in solving Maxwell's equations. When considering symmetries in classical mechanics, 
the discussion is simplified by passing to the Lagrangian or Hamiltonian picture. We present 
a formalism which similarly makes symmetries manifest for Maxwell's equations. In order to 
motivate and illuminate the development of this formalism, we consider the example of light 
propagation in chiral nematic liquid crystals. We believe this to be the first application of Lie 
group techniques to such a problem. The tools we develop, however, are more widely applicable 
to media with a continuous symmetry group. They may also be considered a generalisation 
of 'transformation optics', extending those ideas to allow for the possibility that the dielectric 
tensor and magnetic susceptibility differ. 

We begin in section 2 with a brief discussion of the helical ground state for a chiral nematic 
liquid crystal. In section 3 we introduce the geometry of the optical metric of Joets and Ribotta 
[3], describing the propagation of the extraordinary ray in a uniaxial birefringent material. In 
particular we study this metric and its geodesies for the helical ground state of a chiral nematic. 
The resulting metric is invariant under the simply transitive action of three-dimensional group of 
isometrics which is locally isometric to the Euclidean group E(2) of the plane, which we discuss 
in detail in section 3.1. In fact the isometry is the universal cover E{2) and the Lie algebra is 
of type VIIo in Bianchi's classification. We extend the discussion to Vllt in an appendix. The 
identification of the symmetry group for the chiral phase permits a fully geometrical discussion 
of electromagnetic phenomena, an approach we exploit. The high degree of symmetry allows us 
to solve the Hamilton- Jacobi and wave equations up to one quadrature in the former and up to 
solutions of Mathieu's equation in the latter case. This opens up the possibility of a detailed 
analytic investigation of the type of caustics and optical singularities which should be observable 
in such systems [4, 5, 6, 7]. We then go beyond the geometric optics approximation to consider 
the full Maxwell equations in Section 4. For a uniaxial material whose director field takes on a 
helical configuration we show how the theory of Lie groups leads to separation of variables for 
these equations. The resulting equations take the form of a generalised Mathieu equation. This 
result is similar to others in the literature [8, 9], but the derivation is fully motivated from the 
inherent symmetries of the problem. 

Throughout the paper, we shall make use of the machinery of differential geometry, in partic- 
ular tangent vectors, differential forms and the Lie derivative. References [1, 2] provide readable 
accounts of these concepts. 



2 Chiral Nematics and their Helical Ground State 

A nematic liquid crystal has an order parameter given by a director or direction field specified 
by a unit vector n = (m, 112, TI3), defined up to a sign n ~ — n, with n • n = 1. A chiral nematic 
has a built in twist, specified by a parameter q, which can be realised as a torsion, which alters 
the usual derivative operator acting on a vector by 

= V » n j + ae i3kn k • (1) 

The Frank-Oseen free energy functional in the one-constant approximation is equivalent, up to 
a boundary term, to: 

F[n] = l -J (|V 9 n| 2 - A(n • n - l))d 3 x , (2) 

where we have added a Lagrange multiplier field A to enforce the constraint that n • n = 1. The 
free energy would be minimised if 

V?nj = 0. (3) 
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However, as can be seen by taking another V 3 derivative and skew symmetrising, this is not 
possible over an extended region so the system is frustrated and must adopt some compromise 
configuration [10, 11, 12]. 

One such configuration is the helical phase for which 

n = i cos(pz) + j sin(pz) . (4) 

For more details about the liquid crystals the reader may consult [13, 14, 15, 16, 17, 18]. For the 
helical phase, 

V-n = 0, Vxn = -pn, V 2 n=-p 2 n, (5) 

This configuration is a stationary point of the free energy, satisfying the second order Euler- 
Lagrange equation resulting from extremising F[n]: 

-V 2 n+2gVxn=(A-2 9 2 )n. (6) 

provided we choose 

\=p 2 -2pq + 2q 2 . (7) 

Among these solutions of the Euler-Lagrange equations, the one minimising the free energy 
density has p = q, A = q 2 . The only non- vanishing components of Vfrij are 

V|n 3 = qni Vfn 3 = -qn 2 . (8) 

It follows that the helical ground state is not a solution of the first order frustrated "Bogomolnyi 
equation" 

VS = 0. (9) 

Another means of relieving the frustration is the 'double twist' structure, given in cylindrical 
polar coordinates (p, z, (j>) by 

n = e z cos qp — sin qp. (10) 

Along the z axis, this configuration has = 0, so inside a sufficiently small cylinder, the 

free energy density is in fact lower than that for the helical phase. A structure composed 
of these tubes can fill space, but there will necessarily be defects where the tubes meet [12]. 
Such a configuration gives the so called 'blue phase'. Whether the blue or the helical phase is 
thcrmodynamically preferred depends on the energetic cost associated to accommodating the 
defects of the blue phase. 



3 Optical metrics for Nematics 

If n is the director, and t = ^ , with ds 2 = dx 2 is the unit tangent vector, then then the inverse 
speed or slowness of an extraordinary ray is is given by [3] 

n= ^n2(t.n) 2 +n2(t-n(n-t)) 2 , (11) 

where n a is the refractive index of the ordinary ray and n e that of the extra- ordinary ray. Fermat's 
principle reads 

5 J nds = 0. (12) 
Thus the rays are geodesies of the Joets-Ribotta metric 

ds 2 = n 2 e dx 2 + (n 2 - n 2 e )(n.dx) 2 . (13) 
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Assuming the refractive indices are constants, we can write down the metric for the helical 
ground state given above 

ds 2 — (n 2 cos 2 (pz) + n 2 sin 2 (pz))dx 2 + (n 2 sin 2 (pz) + n 2 cos 2 (pz)) dy 2 

+(n 2 - n 2 ) sin(2pz)dxdy + n\dz 2 (14) 

As another example, consider a particular case where the director field, n, is in a 'hedgehog' 
configuration (cf. [19, 20, 21]) and where in addition the refractive indices n ,n e , vary with 
position inside the ball r = |x| < 1 according to 

n=-, «e = ^=, n °=T^- (15) 
The resulting Joets-Ribotta metric is 

, 9 dx 2 (x.cZx) 2 , . 

which one recognises as that of Hyperbolic Space H 3 in Beltrami coordinates. Remarkably, 
because Hyperbolic space is projectively flat in these coordinates, the light rays are straight lines 
in this case. One could consider the 3-sphere S 3 , by changing the sign in front of r 2 . Then one 
has a model related to Maxwell's fish-eye lens by a coordinate transformation, but whose rays 
are straight lines. Other examples may be found in [19, 20, 21]. Optical properties such as (15) 
my appear unnatural, but modern meta-materials are increasingly able to mimic such refractive 
indices, at least within a certain range of frequencies for the electromagnetic field. 

3.1 E(3) and Left-invariant metrics 

The aim of the present section is to obtain the isometry group of the apparently quite complicated 
metric (14). An enormous simplification results if we use the formalism of metrics on Lie groups. 
We will start with a brief discussion, tailored to the Euclidean group E{2) of isometries of 
the plane, consisting of translations and rotations in two dimensions. Those familiar with the 
construction of left- and right-invariant forms on Lie groups may wish to skip to the summary 
at the end of this subsection. 

We can realise elements of E{2) as matrices as follows. We first fix a number p. To any point 
(X, Y) in the plane, we associate the column vector 



{X,Y)~ \ Y . (17) 




We first note that a general isometry of the plane can be decomposed into a clockwise rotation 
of angle pz about the origin, followed by a translation of (x,y). We can write down a matrix 
M(x,y,z) depending on parameters (x,y,z) which performs this operation as follows 



= M(x,y,z) \ Y . (18) 






Corresponding to any three numbers (x,y, z), we have a unique isometry and conversely each 
isometry corresponds to a unique (x, y, z), provided that we regard z and z + 2ir/p as the same 1 . 
We refer to (x, y, z) as coordinates on the group E{2). 



1 for most of the rest of the paper, we shall in fact work with the covering group E(2) obtained by dropping 
the identification of the z coordinate. This is a slightly technical point which we shall not labour. 
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Our aim is to construct a set of vector fields which are invariant under an action of the group 
E{2). We'll consider for the moment a matrix Lie group G, i-e. a group whose elements are n x n 
matrices for some appropriate n where the group action is by matrix multiplication 2 . We note 
that any element A of Q gives rise to two natural transformations on the group itself. Acting on 
M, a general element of G, they give: 

A A (M) = AM, P A {M) = MA (19) 

and are known respectively as left and right translation, or the left and right action of A. We 
note that the left action and the right action commute: 

A A (P B (M)) = A A (MB) = AMB = (A A (M))B = P B (A A (M)). (20) 

Now suppose that A is infinitesimally close to the identity matrix, A = I+eSA for an infinitesimal 
parameter e. We have 

A A {M) = (I + eSA)M = M + e5 L M (21) 

where the infinitesimal generator S\M = (SA)M can be interpreted as a tangent vector to 
Q at the point M, where we think of Q as a submanifold, i.e. a surface 3 in the space of all 
n x n matrices. This is a right-invariant vector field, since S\(Pb(M)) = Pb(S\M) because of 
equation (20) above. Following a similar procedure, we find that the vector fields SpM = M(5A) 
are left-invariant. 

In order to construct the left- and right-invariant vector fields of E(2), we must first find all 
suitable d A such that / + eSA is an element of E(2). We can do this by using the coordinate 
representation M(x, y, z) above. Since M(0, 0, 0) = /, we can find the most general 5 A by Taylor 
expanding M(x,y,z) for small (x,y,z). We find that the general 6 A is a linear combination of 
the matrices: 
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Mi = M 2 = 1 M 3 = 1 . (22) 



We'd like to express the vector fields in terms of the (x, y, z) coordinates. To do this, we define 
a basis of vector fields for E{2) as follows 

= dM{x,y,z) d = dM(x,y,z) d = dM(x,y,z) 
dx dx ' dy dy dz dz 

We can either think of d/dx as a concrete matrix tangent to E(2) as a surface in the space of 
3x3 matrices, or more abstractly as the vector field which generates a shift from (x, y, z) to 
(x + Sx, y, z) in the coordinate space. This notation captures the fact that under a change of 
variables for the coordinate space, (x, y, z) — > (x 1 , y' , z'), the vector fields transform in the same 
way as differential operators following the chain rule. We can readily calculate the left invariant 
vector field corresponding to Mi : 



cos pz 
shi£ 




L 1 = MM 1 = I sinpz I = cospz— + sinpz — . (24) 
' ' ox ay 



2 Therc are some further assumptions on the smoothness of the group, but we'll take these as given. 
3 We do not assume that a surface is necessarily 2-dimensions, merely that it is of lower dimension than the 
space in which it lives. 
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In a similar way, we can find the rest of the left- and right-invariant vector fields, L, = MM;, Ri — 
MiM: 

L x = cospz-^+smpz^, Rx = 

L 2 = cospzJ^ - sinpz^, R 2 = ^, (25) 

j 1 d n 1 _9 I ™J5 , 9 

^3 — pdz , 3 p dz dy & dx' 

Now that we have the left- and right-invariant vector fields, we can construct the left- and 
right-invariant one-forms which are dual to them. Taking dx, dy, dz to be the one-forms dual to 
d/dx,d/dy,d/dz, the left-invariant forms X 1 and right invariant forms p % are: 



A 1 
A 2 
A 3 



cos(pz)dx + sin(pz)dy, 
cos(pz)dy — s\xi(pz)dx, 
pdz, 



dx + pydz, 
dy — pxdz 
pdz. 



(26) 



The matrices Mi have a natural commutator algebra, the Lie algebra e(2). This determines 
in a natural way the Lie algebra of the vector fields Li and Ri : 



La,, Li 
L3, L 2 



= 0, 
= +L 2 , 
= ~Li, 



and the Maurer-Cartan algebra of the one-forms 



dX 1 
dX 2 
dX 3 



+A 3 A A 2 . 
-A 3 A A 1 . 
0, 



Ri,R 2 

i? 3 , i?i 
R3, R2 



dp 1 
dp 2 
dp 3 



0, 

—R2, 
+Ri, 



-p 3 A p 2 
+p 3 A p 1 
0. 



(27) 



(28) 



The fact that left and right actions commute (20) is reflected in the fact that [L i} Rj} = 0. 
We can check the claimed invariance explicitly. First we note the matrix identity 

M(x, y, z)M{y, i], () = M(x + vcospz — -qsmpz, y + ijcospz + vsinpz, z + Q, (29) 

from which we deduce that the element M{v,r],(,) acting by right translation takes (x,y,z) to 
(x', y',z') where 



x = x + vcospz — r/sinpz, 
y' = y + -qcospz + vs'mpz, 

z' = z + C 

Making these substitutions into p 1 , treating i], ( as constants, we find 

dx+pydz = dx'+py'dz', 
dy — pxdz = dy'—px'dz', 
pdz = pdz', 



(30) 
(31) 
(32) 



(33) 
(34) 
(35) 



so that the p % are indeed invariant under right translations. Interchanging the roles of x, y, z and 
v, r\, £ in (29) we deduce that the element M{v,rj,C) acting by left translation takes (x,y,z) to 
(x', y',z') where 



= xcosp( — ysinpC + v , 
= ycosp( + xsinp( + i], 
= z + C- 



(36) 
(37) 
(38) 



G 



substituting into X 1 , again treating v, rj, £ as constants, we can check that 



cos(pz)dx + sm(pz)dy = cos(pz')dx' + sm{pz')dy' , (39) 
cos(pz)dy — ain(pz)dx = coa(pz / )dy / — sin(pz')c?a;', (40) 
pdz = pdz! ', (41) 

so the \ l are invariant under left translations. 

Armed with these invariant one-forms, we are now in a position to construct metrics which 
are invariant under an action of E(2). For example, the flat metric can be written in terms of 
the left-invariant one-forms as: 

d s 2 =p- 2 (A 3 ) 2 + (A 1 ) 2 + (A 2 ) 2 , (42) 

and is hence manifestly left-invariant. Now for the helical ground state of the nematic liquid 
crystal, 

n • dx = A 1 , (43) 
thus the Joets-Ribotta metric of the helical phase may be written as 

dsl = n 2 (p- 2 (A 3 ) 2 + (A 1 ) 2 + (A 2 ) 2 )+(n 2 -n 2 )(A 1 ) 2 , 

= n 2 (A 1 ) 2 + n 2 (A 2 ) 2 + ^f(A 3 ) 2 (44) 

which is a left invariant metric on E(2). In fact, any left invariant metric may be brought into 
this form by a global right action of E(2). Cartan's formula for a p-form reads 

C x w = ixdu + d{ixw) , (45) 

and so the non-vanishing Lie derivatives are 

£l 3 A 1 = -A 2 , (46) 

C L3 \ 2 = A 1 , (47) 

C Ll X 2 = -A 3 , (48) 

Cl.X 1 = A 3 . (49) 

Thus while L% is an additional symmetry of the flat metric none of the Li are symmetries of the 
Joets-Ribotta metric. 

To summarise then, (x,y,z) £ M 3 may be considered as coordinates on E(2), the universal 
cover of the two-dimensional Euclidean group E{2) with A 1 left-invariant one-forms and Li left- 
invariant vector fields. If we were to identify the coordinate modulo n = 1, 2, . . . the group 
would be the n -fold cover of the Euclidean group E(2), which corresponds to n = 1. With 
this identification, the Joets-Ribotta metric of the helical phase is a left-invariant metric. Its 
symmetry algebra, the Lie algebra e(2), is of Bianchi type VII . As an aside, this may be obtained 
from the rotation group algebra so (3) by means of a Wigner-Inonii contraction. If Mj are the 
generators of so (3), one sets Mi = \M\ , M 2 = \M 2 , M 3 = M 3 and takes the limit e -> 0. Under 
this contraction, the direction field n-cix = L\ arises as the image of the Hopf fibration generated 
by the right action of L\ and the Joets-Ribotta metric as the image of the Berger-Sphere. For 
the relevance of the Hopf fibration to chiral nematics see [10, 11, 12, 16, 22, 23, 24, 25]. 
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Figure 1: Two geodesies of the Joets-Ribotta metric (14). The curves are shown projected into 
the x — z plane. One is unbounded in z, whereas the other is bounded. Not shown is the motion 
in the y direction which gives both these curves a cork-screw motion. 



3.2 The Ray Approximation 

In the ray approximation we are looking at geodesies with respect to a left-invariant metric on the 
universal cover of the Euclidean group E(2). Now the Euclidean group E(2) is the configuration 
space for a rigid body in two-dimensional Euclidean space E 2 . The motion of a rigid body moving 
in a homogeneous, incompressible, inviscid fluid [26] is known to correspond to geodesic motion 
with respect to a left-invariant metric on the Euclidean group [27, 28]. The present situation 
corresponds to a cylinder with its axis in a plane [26, 29] which may be reduced to the quadrantal 
pendulum. 

To see this in detail note that the Eikonal equation is 



§( L 3^) 2 + ^(W) 2 + ^(iMr;: ,A (r,(» 



and it separates. If W = k x x + k y y + G(z), then 



-o(-r-f + \( k x cos(pz) + k y sin(pz)) 2 + \(k x sm(pz) - k y cos(pz)) 2 = u 2 . (51) 
ni az nt n z 



The Killing vectors Ri give rise to three constants of the motion of the form 

Pi =9^1%. (52) 

of which two, pi and p 2 mutually commute. 

We may immediately find the equations for rays in a first order form by making use of the 
relation: 

^-=g^^L. (53) 

dt y dx» K ' 



We find: 



<± - ^ U - ^)<-(2,k -r). (54) 
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nli - + (55) 

..^ = cX-y(l+||) + y(l-|)cos(2pz-«). (56) 

Here tan?/> = k y /k x , k — ^Jk 2 + k 2 and tan# = 2k x k y /(k 2 — k 2 ). Let us first consider (56). 
Introducing new constants a, (3 and defining ( = pz — 9/2, we find 

C 2 -4(" + /W2C))=0 (57) 

which is the so-called quadrantal pendulum equation [29] . The pendulum has two different types 
of behaviour, depending on the constants a and /3. If a > |/3|, then |£| and hence \z\ will increase 
without bound. This corresponds to a pendulum swinging through complete revolutions. If 
a < \j3\, then £ will oscillate about 2n-n for some integer n. This corresponds to the standard 
libratory motion of a pendulum. Thus we find two behaviours for the rays. Either the rays can 
penetrate in the z-direction or else they are trapped to move between two planes perpendicular 
to the z-axis. Finally, we can consider the other equations of motion, (54) and (55). We can 
interpret these as saying that the tangent vector of the ray oscillates around an average direction. 
For rays which are not bounded in z, the result is a 'cork-screw' curve, similar to a helix. The 
'tightness' of the spiral is determined by how close n 2 /n 2 is to 1. Figure 1 shows some examples. 



3.3 The wave equation 

For an approximate description of light propagation beyond the ray approximation one may 
use the scalar wave equation. The scalar wave equation captures some of the wave aspects of 
light, whilst ignoring the complications relating to polarisation which arise in the full Maxwell 
equations. We may expect the wave equation on the Joets-Ribotta metric to share some features 
with one of the two polarisations of the full Maxwell equations. It takes the form: 

d 2 ^ n 2 1 1 

= -^t Y + ^L 3 L 3 <!> + — L 2 L 2 ^+— LiLi*. (58) 

0% TLa Tip 



(oj 2 n(k x cos(pz) + k v sin(pz)) 2 n(kx sin(pz) — k v cos(pz)) 2 \ F = . (59) 

V ni ni J 



It separates. That is if * = e i( - +k - x+k y y - w ^F(z), then 

1 d 2 F ( o 1 o 1 

n 2 dz 2 v rig 

Recalling a, f3, £ from the previous section, this is of the form 

^ + (a + /3cos(2C))F = 0, (60) 

which is Mathieu's equation. 

By the Floquct-Bloch theorem, the general solution of (60) is of the form 

F = Cl e^/(C) + c 2 e-^/(-C) (61) 

where f(() = f(( + 2n) and p depends on a and (3. Expanding /(£) as a Fourier series, we 
deduce the Laue-Bragg conditions that an incoming wave with wave vector kj n incident on some 
region where propagation is described by (58) is reflected/diffracted with wave vector k out where 

p(k out - k in ) z = meZ. (62) 
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We may think of \x = fj,(k x ,k y ,u>) as defining a dispersion relation, averaged over the period 
in the vertical direction. When /x is real we expect propagating waves, whereas when /j has an 
imaginary component the solutions either decay or grow exponentially in £. It can be shown that 
the marginal cases between propagation and damping occur when /i = 0, tt (note /i is only defined 
up to multiples of 2ir). This defines a set of surfaces in the (k x ,k y ,uj) space which separate out 
the regions where the wave propagates and where it is damped. To determine these surfaces, we 
can (for \x = 0, the other case follows similarly) expand F in Fourier series: 

oo 

F = ^c„e^ (63) 

— oo 

and obtain a three term recurrence relation 

~n 2 c„ + ac n + 2f3(c n - 2 + c n+2 ) = . (64) 

The condition that this relation admits a non-trivial solution may be related to the vanishing of 
an infinite determinant, a procedure known as Hill's method, see e.g. [30]. 



4 Maxwell's equations 

Before we discuss Maxwell's equations for the helical phase of a nematic liquid crystal, we will 
first formulate Maxwell's equations for a general medium in the language of differential forms. 
This will be the most convenient language in which to discuss how to apply the machinery of 
Lie groups to the problem in hand. We will work directly with the fields rather than introducing 
potentials as this avoids tackling the issue of gauge invariance. We'll work on a 4-dimensional 
manifold M, but this restriction is not necessary. We begin by noting that we can define two 
2- forms: 



= ^F^dx* A dx v = E i dx t Adt+ ^e ilk B l dx° ' A dx k , 
1 „ , „..-.,. 1 



G = -G^dxP A dx v = H t dx l Adt- -t iik D l dx 3 A dx k . 

These forms encode the fields contained in the antisymmetric 4-tensors with components 

Fio = Ei Fij — eijkB k 

Gm = Hi Gij = —Eij^D 

Where Ei, B % are the electric field and magnetic displacement, D l is the electric displacement 
field and Hi is the magnetic field. The advantage of packaging the fields as two-forms is that 
Maxwell's equations become the simple pair of relations 

dF = 0, dG = J, (66) 

with J the current 3-form. Of course, to close this system of equations for a prescribed J, we 
must specify a relation between F and G, the constitutive relation. This is nothing more than 
the usual relations one requires relating (E, B) and (D, H). In the language of forms, we require 
a map from the space of sections of VL 2 (M) to itself 4 . In many materials, the constitutive relation 



4 In higher dimensions, G will be a n — 2 form, but similar considerations apply 
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is local and linear, so may be represented by a section of the bundle End(Q 2 (M), i.e. a possibly 
space dependent linear map C such that 



G = CF, i.e. G^ = C^ T F KT (67) 

where C acts pointwise. This tensor C, together with the differentiable structure of the manifold, 
is the minimal data required to define Maxwell's equations - it has not thus far been necessary 
to introduce a metric or other structure to M. In index notation the Maxwell equations take the 
form: 

dfrF va] = 0, 8 [IX (C„ a] KT F KT ) = 0. (68) 

We note that defining C as an endomorphism, i.e. with two indices up and two down, ensures 
that it is not necessary to define a connection in order to take derivatives covariantly. 

In order that (66) defines a suitable hyperbolic system of partial differential equations, re- 
strictions are required on C. We shall assume that C satisfies some such suitable conditions, 
without specifying what those might be. As a simple example, we may take C to be the Hodge 
map induced by a Lorentzian metric g, i.e. we take C — * g . If g is the flat metric, this gives the 
classical Maxwell equations in the vacuum. If g is not flat, we may interpret the field as an elec- 
tromagnetic field propagating in a gravitational background. In the case that g is static, we can 
alternatively interpret the field as propagating through some material with a position dependent 
dielectric tensor, and magnetic susceptibility fiij. This is the basis of transformation optics 
[31, 32, 33]. For a material in which Maxwell's equations have a gravitational interpretation it 
must be the case that = /i^ in suitable units, i.e. the material is impedance matched [33]. 
This need not be the case for a general material. We will take C to have the following form: 

C(dx l A dt) = - ]^Si ej kl dx k A dx l 
C(dx i A dx j ) = e^ k {^ x ) k idx l A dt. (69) 

Note that if C = * s for some Lorentzian metric, we have C 2 = —1 so that £ij(/J.~ 1 )jk — Sj k , 
justifying our assertion that materials with a gravitational analogue are impedance matched. 

The liquid crystals in which we are interested are uniaxial, so that at each point, we can 
assume that the quadrics defined by e and fj, are spheroidal (i.e. ellipsoids with an axis of sym- 
metry) with a common axis. In other words, there is locally a basis in which the tensors have 
the form 

' s\\ 

e ± ) n = | fj, ± | . (70) 


If we assume that the axis of the material lies along n, this can be written in a more covariant 
form as 

£ij = e±Sij + (e|| - e±)riinj, = fj,±8ij + (fi\\ - fj,±)riinj. (71) 

Before we discuss the consequences of such a constitutive relation in the case of a nematic 
liquid crystal in the helical ground state, let us first consider for a moment the geometric optics 
approximation. 

4.1 Geometric optics 

Let us consider the Maxwell equations described above in a geometric optics limit. We consider 
a field which takes the form 

F = e^ (F + aF 1 + ...) (72) 



Mil 

























11 



where by assumption Fj are 0(1) as a — > 0. We assume that there are no currents or charges, 
so that Maxwell's equations become dF = dG = 0. We also assume that C varies slowly by 
comparison to the wavelength of the field. Inserting our ansatz and collecting terms in a, we 
find 



oo 

I 



= -dS AF + J2( dF k-i+idS /\F k )a k -\ 



a 

k=l 

oo 

= -dS ACF + J2( dCF k-i+idS ACF k )a k -\ (73) 

a k=l 

Let us first consider the 0(a~ 1 ) terms. This is a system which asserts that F is in the kernel of 
a linear operator which maps from one 6-dimensional space to another 6-dimensional space. The 
condition that a non-trivial Fq exists gives a differential condition on S involving C which we 
interpret as the eikonal equation. Associated to a solution of the eikonal equation is a 2-form F 
which gives the polarisation of the wave. In general, there will be only one polarisation associated 
to each solution of the eikonal equation. Once we have solved for S and Fo, we can inductively 
construct F k by solving the equations 

= dF k -i + idS A F k , 

= dCFh-i + idS A CF k . (74) 

Presumably the well posedness of this system is a necessary condition that C be an acceptable 
constitutive map. 

In the case where C = * g , the eikonal equation can be shown to reduce to 

dS A * g dS = 0, (75) 

which is the Hamilton-Jacobi equation for geodesies of the metric. In this case, there is a two 
dimensional space of possible polarisation tensors. They take the form 

F =dSAfo, g(dSJ o ) = 0. (76) 

The Hamilton-Jacobi equation requires that dS be null. Suppose for example that at a point, 
dS is parallel to dt — dx, then the space of polarisations at that point is spanned by dS A dy and 
dS A dz. 

In the case where C has the uniaxial form introduced above, the eikonal equation reduces to 
the form: 

-^S 2 t + -VS 2 + ( — - -)(n ■ VS) 2 ) (-e ± S 2 + —VS 2 +(±-±)( n - VS) 2 ) = 

(77) 

The medium is thus birefringent. We see straight away that the condition on S factors into two 
separate Hamilton-Jacobi equations associated to the two metrics 

dt 2 

9b = +£||cZx 2 + (£_ L -£||)(n-dx) 2 , (78) 

dt 2 

g E = h ^||dx 2 + (fj, ± - M||)(n • dx) 2 . (79) 

These are both of the Joets-Ribotta form we have previously considered. It can be checked 
that the polarisation tensor associated to a solution of the Hamilton-Jacobi equation of gs 
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has eijkFijUk = B n = 0, whereas for a solution of the Hamilton-Jacobi equation of qe, the 
polarisation tensor has riiF it = E n = 0. Note that we do not require that n remains constant 
for this derivation, provided it varies slowly compared to the wavelength of the light. In the case 
that n varies from point to point, the polarisation will also change so that to leading order in 
a, either the magnetic or electric field parallel to the director will vanish, depending on which 
type of ray we consider. Often, one takes = fJ*\\ in which case, gE is simply the Minkowski 
metric and its geodesies are the ordinary rays. The rays of the metric gs are the extraordinary 
rays and gs is the Joets-Ribotta metric, where we identify e±fi± = n 2 Q and £||/Uj_ = n\. 

If C is of the form (69), but with no uniaxial assumption, then the rays will typically be 
geodesies of a Finsler geometry. 

4.2 Symmetries 

So far, we have re-cast familiar results into the notation of differential forms. Whilst this is a 
satisfying exercise, it is not clear that it introduces any benefits beyond putting the equations 
in a manifestly coordinate invariant form. For our purposes, the great advantage is that this 
form of the equations permits a concise discussion of the symmetries of the system and allows 
the machinery Lie groups to be brought to bear. We start by defining a Killing vector K to be 
a vector which satisfies 

C K C = 0. (80) 

Recall that C is simply a tensor, so the Lie derivative is defined as a consequence of the diffcr- 
cntiable structure of M. Making use of this and Cartan's relation, we deduce that if a 2-form F 
obeys Maxwell's equations: 

dF = 0, d(CF) = 0, (81) 

then so will CkF and in particular, the diffcomorphism induced by K will map solutions of the 
equations into solutions of the equations. An important example occurs when C = * g and K is 
a Killing vector of g. 

Suppose that we have a group which acts simply transitively on M by left actions and which 
preserves the material configuration, as is the case for the E{2) x R t symmetry of helical ground 
state of the nematic liquid crystal. Then it must be that C may be written in terms of the left 
invariant one-forms and their duals as: 

C= l -C ab cd {\ a h\ b )® (L c AL d ) (82) 

where C a b cd are some constant coefficients. Here, indices run over 0, . . . , 3. We can make use of 
this to write down Maxwell's equations for a nematic liquid crystal in its helical state. We take 

F = E i X i Adt+^kB^ AA fe . (83) 

This choice of basis is very similar to the rotating basis chosen by Peterson, who investigated 
the electromagnetic field propagating through a nematic liquid crystal in its ground state [8] . In 
our case, this choice of basis arises naturally from the group structure of underlying symmetries. 
We assume further that 

C(\ l Adt) = - l - £lJ e okl \ k AX 1 
C(A*AA J ') = e ljk ( f i- 1 ) kl X l Adt. (84) 
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Where e, [i have the uniaxial form we previously assumed (70). Maxwell's equations for the 
electric and magnetic fields take the form 



e ijk Lj(E k ) + - PijEj - 0, 



at 



P%j{^ 1 )jkBk — J i- 



(85) 



(86) 



Li(Bi) = 0, 

£ijLi(Ej) = p, (fj,~ 1 ) k ie i jhL j (B l ) - e iy 

The matrix P^j has non-zero components 

Pn = P22 = 1- 
These equations can be separated with the ansatz 

E . = e i(k x x+k y y- U t) f.( z ^ B . = e i(k x x+k y y-wt) g ^ z y ( 87 ) 

The components fz{z), gz{z) are given by a linear combination of other components, so that the 
Maxwell equations reduce to a system of differential equations of the form: 

F'(z) + (a + p ie 2ipz + l3 2 e- 2lpz )F(z) = 0. (88) 

Here F(z) = (fi(z), f 2 (z), gi(z), 32(2))* is a 4- vector and a, fa are 4x4 matrices, given by: 



a = 





-1 


ip|i|fc| 2 . 



»N 2 

2bj 



1 






- 4 |k| 2 



»N 2 



2£±fl±LO 





— ioj 







(89) 



and 



K 

4cj 



/ 


1 






i 

.til 



— i 











(90) 



where we have introduced n = k x + ik y . We see that the Euclidean symmetry of the original 
problem is still manifest since a rotation in the x-y plane sends k — > e l6 K, which is cancelled 
by a suitable shift in the z coordinate. We may view (88) as a generalised Mathieu equation. 
Mathieu's equation itself may be written in this form with 2x2 matrices. By Floquet's theorem, 
the general solution of (88) will take the form: 



F(z) 



z h 2 {z) 



z h 3 (z) 



z h 4 {z) 



(91) 



where hi(z) = hi(z + ir/p) are 4- vectors. Making use of discrete symmetries of the equations, one 
may show that if ji is a Floquet exponent, then so is — \i and ~p, implying relations amongst the 
fj,i. This equation may be studied using the infinite determinant techniques of Hill, an approach 
similar to that of [9], but that takes us beyond the scope of the current paper. We hope to address 
this issue in future work. Since we have retained the independence of the magnetic susceptibility 
and the permittivity, this analysis applies equally well to magnetic materials with helical phases 
[34]. 
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5 Conclusion 



We have shown that certain properties of the chiral phase of a nematics liquid crystal are inti- 
mately tied to the symmetries it possesses. We have shown that the Joets-Ribotta metric, which 
describes the propagation of extraordinary rays, is a left-invariant metric on E(2) and we have 
shown how the underlying symmetry group can be practically used to understand properties of 
waves in such a medium. 

We have separated the Hamilton- Jacobi equation and the wave equation for this metric. The 
wave equation can be reduced to Mathieu's equation and the Hamilton- Jacobi equation to the 
quadrantal pendulum equation. We have also seen how Maxwell's equations for a general uniaxial 
material whose director field lies in a helical configuration can be reduced to coupled ordinary 
differential equations generalising Mathieu's equation via a novel application of the theory of 
Lie groups. This new formalism is applicable to the macroscopic Maxwell equations whenever 
the medium has a continuous symmetry group. The approach taken generalises transformation 
optics to permit non-impedance matched media. 

As we have seen even in this simple example, the extraordinary light rays propagating through 
a liquid crystal explore a much richer geometry than the usual flat geometry of light rays in the 
vacuum. This opens up the possibility of constructing analogues for the propagation of light in 
a gravitational field. In this case the light rays in the liquid crystal may be mapped onto light 
rays propagating in a Bianchi VIIq cosmology [35] whose spatial sections have a fixed geometry, 
but one may imagine more ambitious possibilities. 



A Generalization to Bianchi type Vllh 

It is interesting to ask whether the set up above generalises to the Bianchi type Vllh group. For 
this section we set p = 1, in order not to clutter up the formulae. 
We now define left-invariant one-forms and dual vector fields by 

A 3 = dz, L ^d~z' (92) 

d d 

A 1 = e hz (cos zdx + smzdy) , L± = e~ hz (cos z— + sin z — ) , (93) 

Ox ay 

d d 

A 2 = e hz (cos zdy - sin zdx) , L 2 = e (cos z- sinz)— . (94) 

ay ox 

The right-invariant one-forms and vectors fields are 

p3 = dz > R * = Wz +x Wy- y Tx- h{x o\ +y ib' (95) 

d 

p 1 = dx + (1 + h)ydz , i?i = ^ , (96) 

d 

p 2 = dy-(l- h)xdz , R 2 = — . (97) 



dy 



The metric 5 



(A 2 + A 2 + A 2 ) = n\ (dz 2 + e 2hz (dx 2 + dy 2 )) (98) 



3 In what follows n e and n will be taken to be constant, that is position independent. 
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is in fact that of hyperbolic three space, in the upper half space or Poincare patch space model. 
Setting 

e* = ± y=\ (99) 

it becomes 

2 

^(dz 2 + dx 2 + dr 2 ) , (ioo) 

and we see that optically we can think of a vertically stratified isotropic medium with Cartesian 
coordinates (X, Y, Z) and refractive index 

E- (101) 

Rays are now circles orthogonal to the plane Z = 0. 
The metric 

ds 2 = n\ (A 2 + A 2 + A 2 ) + (n 2 - n 2 )A 2 (102) 

may thought of as describing a vertical stratified anisotropic medium with extraordinary and 
ordinary refractive indices varying with height Z in the same way, i.e. as 

, and respectively . (103) 

Such a variation might be due to temperature variation within the material, for example. As 
before, the wave equation separates but F(z) now satisfies 

d 2 F dF ' -2hz m 2 a -2hz 



dz 2 



+2h — +( w 2 n 2 -^(l + ^|)(A : 2 +A : 2 ) + ^(fc 2 +fc 2 )(l-^|)cos(2z-0))F^O. (104) 
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